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ABSTRACT 

We predict the formation histories, properties and descendants of Lyman-break galaxies 
(LBGs) in the ACDM cosmology. In our model, which incorporates a top-heavy IMF in star- 
bursts, we find that most LBGs are starbursts triggered by minor mergers of galaxies. The 
duration of the LEG phase is predicted to be quite short, ^ 20 — 60 Myr We investigate the 
distributions of stellar and halo masses and morphologies for bright (Luv > L^y) and faint 



uv> 

(Luv > O.lL^y) LBGs at z = Z, z — Q and z — 10 (where we classify LBGs according 
to their rest-frame UV luminosities relative the observed characteristic luminosity L'^y at 
z ~ 3). Bright LBGs at z = 3 are predicted to have median stellar masses ^ 1 x l{fh^^MQ 
and host halo masses ^ 3 x IQ^^h^^ Mq, and to be typically mildly disk-dominated in stellar 
mass. On the other hand, faint LBGs at z = 10 are predicted to have median stellar masses of 
1 X l{)^h^^MQ and host halo masses ~ 2 x lO^^ft-^^A/©, and to be generally bulge- 



only ' 

dominated. Bright LBGs at z = 3 evolve to present-day galaxies with median stellar mass 
^ 5 X IQ^'^h^^MQ (comparable to the Milky Way), consisting of roughly equal numbers of 
disk- and bulge-dominated systems, and hosted by halos with median mass ^ 2 x IQ^^h^^MQ 
(corresponding to medium-size galaxy groups). The model predicts that 40% of Milky Way 
mass galaxies at the present-day have a bright LBG progenitor in the redshift range 3 < z < 4, 
while 95% have a faint LBG progenitor in the same redshift range, and 7% have a faint LBG 
progenitor atlO<z<ll. With our multiwavelength model, we also investigate the overlap 
between the LBG population and that of submillimetre selected galaxies (SMGs); at z = 3, 
only ^ 1% of bright LBGs are also predicted to also be bright SMGs (with an 850/im flux in 
excess of 5 mJy). 
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1 INTRODUCTION 

High redshift galaxy populations identified using the Lyman-break 
drop-out technique have played a pivotal role in the development 
of our un derstanding of the global star formati on history of the 
Universe jMadau et al.1ll996l : isteidel et al.lll999| ). The star forma- 
tion rate density inferred from measurements of the galaxy lumi- 
nosity function in the rest frame ultraviolet (UV) displays a broad 
peak around redshift ~ 2, p rior to a steep decline towards the 
present day jLillv et al.lll99q) . This type of behaviour, which cor- 
responds to a gradual build up of the global stellar mass in galax- 
ies, arises naturally in h ierar chical structure formation models, as 
shown by ICole et al.l jl994h and confirmed in a ran ge of calcu- 
lations with varying galaxy formation physics (e. g. iBaugh et aTl 



1 19981 : ISomerville et aUl200ll : ISpringel & Hemquislll2003h . In gen 
eral terms, the star formation rate density is set by the abundance 
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of haloes as a function of mass which are in place at a given epoch, 
the ability of gas to cool within these haloes and the effectiveness 
of supemovae at reheating gas which has cooled. The star forma- 
tion density peaks around z ~ 2 because gas can cool efficiently at 
high densities and stays in the cold phase because there is, for the 
first time, a large population of haloes in which supernova heating 
is ineffective. In this paper we investigate the properties of Lyman- 
break galaxies (hereafter LBGs) in a hierarchical galaxy formation 
model and show how they are related to galaxies in the present day 
Universe. 

The Lyman break technique is based on observing a break 
bluewards of the Lyman limit at 912A in the rest-frame spec- 
tral energy distribution of a galaxy. The flux at these wavelengths 
is expected to be small for a number of reasons; i) few pho- 
tons can be produced by stars at these short wavelengths, ii) such 
energetic photons can be absorbed by the atomic hydrogen in 
the interstellar medium within the galaxy, and iii) the photons 
can also be absorbed by intergalactic neutral hydrogen along the 



© 2009 RAS 



2 J. E. Gonzalez et al. 



line of sight from tfie galaxy to the Earth. Galaxies identified 
with t his t echnique are known as Lyman-break galaxies (hereafter 
LBGs^idel & Hamilton 1992; Steidel et al. 1996). The Lyman- 
break technique uses multi-band photometry. By using a suitable 
set of filters, galaxies in a particular range of redshifts can be 
isolated. Ground-based observations are restricted to applying the 
Lyman-break technique for galaxies at z > 2.5. For this redshift 
range, the absorption edge, which occurs in the far-UV, is shifted 
to wavelengths that are not blocked by the Earth's atmosphere. The 
multi-band photometry is also used to make additional colour selec- 
tions to minimize the contamination from foreground objects (both 
stars and lower redshift galaxies). LBGs are selected to be red in 
the shorter wavelength colour, due to the presence of the break, and 
blue in the longer wavelength colour, to show they are star-forming 
galaxies. Hot, high mass stars with short lifetimes tend to dominate 
the emission in the UV. Such stars spend a significant fraction of 
their lifetime within dense giant molecular clouds and their light 
is heavily extincted. Hence in order to infer a galaxy's star forma- 
tion rate from its UV luminosity requires an uncertain and quite 
often large correction for dust extinction. Furthermore, an assump- 
tion about the fom of the initial mass function with which stars are 
produced is also required in order to infer the star formation rate. 

In this paper we use a semi-analytical ap proach to mod el the 
galaxy formation process (seefBaugh 200d and lBensonl20 1 Ol for re- 
views). This type of model can be used to make predictions for the 
properties of high redshift galaxies and can connect these objects 
to their present day descendants. This approach was first applied to 
the Lyman-break population by 'Baugh e t alj ( Il998h . who showed 
that the model of ICole et alj (1994) could rep roduce the observe d 
surface density of Lyma n-break galaxies (e.g. ISteidel et"al]|l996h . 
ISomerville et al.l ( 1200 ll) applied similar techniques to argue that 
bursts triggered by minor and major galaxy mergers played a role 
in producing the LBG population. In this paper we use the version 
of the D urham GALFORM g alaxy formation simulation code pub- 
lished bv Baugh et alj bOOSh to study LBGs and their descendants. 
This model invokes a different form for the stellar initial mass 
function in merger-driven starbursts, in order to reproduce the ob- 
served counts and redshift distribution of galaxies measured at sub- 
millimetre wavelengths. Baugh et al. (2005) showed that this model 
also matches the observed LBG luminosity function at 2 = 3, as 
well a s giving a reasonable m atch to the present day galaxy popu- 
lation. iGonzalez et all j2009l) extended the local tests of this model 
to include the size-luminosity relation for disk-dominated galaxies. 
The semi-analytical model natural ly relates galaxies at h igh red- 
shift to their local counterparts. In iGonzalez et al.l ( 1201 ll) we fol- 
lowed the evolution of submillimetre galaxies (SMGs) and their 
descendants. In this paper, we carry out a similar analysis in the UV 
by studying the evolution and descendan ts of Lyman-brea k galax- 
ies. In a companion paper to this one, L acev et ^ ( I2OIII) present 
model predictions for the evolution of the far-UV luminosity func- 
tions over the redshift range z — 3 — 20, and make a detailed com- 
parison of these with current observational data. They also make 
predictions for the amount of dust extinction, and for the evolution 
of a wide variety of other properties of LBGs over a range of red- 
shifts. The present paper instead looks in more detail at the forma- 
tion histories and lifetimes of LBGs and the distributions of stellar 
and halo masses and morphologies for LBGs at selected redshifts 
and luminosities, and investigates the corresponding properties for 
their present-day descendants. 

This paper is laid out as follows. In Section 2, we outline the 
GALFORM galaxy formation model. In Section 3, we review the 
model predictions for the luminosity function of LBGs, show ex- 



amples of galaxy merger histories which include LBGs, and inves- 
tigate the triggering and duration of the LBG phase. We also present 
the probability that a present day galaxy or halo had an LBG pro- 
genitor. In Section 4, we study different properties of the LBGs, 
such as stellar and host halo masses, and morphologies at redshifts 
z — 3, z = 6 and z = 10, and discuss the dominant mode of star 
formation behind the LBGs. In Section 5, we examine the over- 
lap between the LBG and SMG populations. Finally, we present a 
summary in Section 6. We quote magnitudes on the AB system and 
define the Hubble parameter as h — Hq/WO km s~^ Mpc^^. 



2 THE GALAXY FORMATION MODEL 

We first give a brief overview of the Durham semi-analytical galaxy 
formation code, GALFORM, before discussing two aspects which 
are of particular relevance for this paper, galaxy mergers and the 
form of the IMF adopted in different mo des of star formation. 
Full details of GALFORM can be found in ICole et al.l ilOOd) and 
Benson et al 1 (l2003h . The model used here is described in full in 



Baugh et al]l2005l) and lEacev et alj ( I200&) . 



2.1 Basic Components 

Galaxies are assumed to form inside dark matter haloes, with their 
subsequent evolutio n being controlled by the merger and accretion 
history of the halo dWhite & Rees|[l978h . The physical processes 
modelled in GALFORM include: i) the hierarchical assembly of dark 
matter haloes; ii) the shock heating and virialization of gas inside 
the gravitational potential wells of dark matter haloes; iii) the radia- 
tive cooling of the gas to form a galactic disk; iv) star formation in 
the cool gas; v) the heating and expulsion of cold gas through feed- 
back processes such as stellar winds and supemovae; vi) chemical 
evolution of gas and stars; vii) galaxy mergers within a common 
dark halo as the result of the decay of a galaxy's orbit due to dy- 
namical friction; viii) the evolution of the stellar populations using 
population synthesis models; ix) the extinction and reprocessing of 
starlight by dust. 

We use the s ame galaxy form a tion p arameters adopted by 
iBaugh et al.l bOOSh and iLacev etal] bOOSi) . The model is set in 
a cold dark matter plus cosmological constant cosmology, with 
present-day matter density parameter Q,o = 0.3, a cosmological 
constant Q,\ — 0.7, a baryon density Q,b = 0.04 and a power spec- 
trum normalization given by erg = 0.93. 

In this paper one of our primary goals is to connect high red- 
shift galaxies with their present-day descendants. This requires a 
prescription for building halo merger histories which is accurate 
over a long i nterval in time. To ach ieve this we use the Monte-Carlo 
algorithm o f lParkinson e'taL|(l2008i) . which is a modified version o ' 
that proposed by Cole et al. 1 2000l) and used in lBaughetalJ ( l2005h . 
The IParkinson et al.l recipe has been tuned to mat ch the merger 
histor ies e xtracted from the Millennium Simulation (ISpringel et al.l 
I2OO5I ; see iHarker et al.ll2006 for a description of the construct ion 
of merger trees from the simulation). We note that lBaugh et alj did 
not make any predictions which depended on the accuracy of the 
trees over long time intervals. Rather, they laid down grids of halos 
at a range of redshifts to compute the luminosity function of SMGs 
and other galaxies at these redshifts, and then, by integration, the 
number counts of SMGs a nd the cosmic sta r formation history. We 
have checked that using the lParkinson et aL|tre e algorithm does not 
alter the predictions presented in Baugh et al.l 

Dust extinction is expected to have a large effect on the far-UV 
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luminosities of LBGs. We cal culate dust extinctio n in the model us- 
ing the approach described in* Lacev et al.l ( 1201 lb . and we refer the 
reader to that paper for a discussion of the effects of dust extinc- 
tion on LBGs. In summary, the dust extinction is calculated using 
a two-component dust model, with molecular clouds within which 
stars form embedded in a diffuse dust medium mixed with the stars 
and distributed in an exponential disk. All of the predictions which 
we show in this paper include the effects of dust extinction. 

2.2 Galaxy mergers and merger-driven starbursts 

When dark matter haloes merge, the galaxy in the most massive 
halo is assumed to become the central galaxy in the new halo while 
the other galaxies present become its satellites. In the standard 
form of GALFORM, we assume that central and satellite galaxies 
can form stars from their cold gas reservo irs, but tha t only central 
galaxies can continue to accrete cold gas. jpont et aT 2008 explore 
a variant of this model in which satellites are allowed to retain some 
fraction of their hot halo, depending on their orbit within the main 
halo, with gas allowed to cool onto the satellite disks.) 

In the standard version of the model, mergers are assumed 
to happen only bet ween satellites an d central galaxies within the 
same halo (although iKim et al]|2009 considered a variant of this). 
Bursts of star formation are assumed to be triggered by galaxy 
mergers. Two types of mergers are defined, major mergers and mi- 
nor mergers, according to whether the ratio of mass of the smaller 
to the larger galaxy exceeds the value of the parameter fdup, i.e. 
M2/M1 > /eiiip. Bursts are assumed to be triggered in all major 
mergers, and also in minor mergers which satisfy M2 /Mi > /burst 
and where the ratio of cold gas mass to stellar plus cold gas 
mass in the larger galaxy exceeds /gas.crit . The parameter values 
adopted in the Baugh et. al model are /ciiip = 0.3, /burst=0.05 and 
/gas.crit = 0.75. In a major merger, the stellar disks of the merging 
galaxies are assumed to be transformed into a new stellar spheroid 
(resulting in a pure bulge galaxy) and any gas present in the disks 
is converted into spheroid stars in a burst. In minor mergers, the 
stellar disk of the larger galaxy is preserved and the stars of the 
smaller, accreted galaxy are added to its bulge (resulting in a disk 
plus bulge galaxy). If the condition for a burst in a minor merger is 
satisfied, then the gas of the smaller galaxy is converted into stars. 
If the condition is not satisfied, any accreted gas is added to the gas 
disk of the central galaxy. From this description of the possible sce- 
narios following a galaxy merger, it is readily apparent that mergers 
change the morphology of galaxies in the model, driving the mor- 
phology towards early types (bulge dominated), with the cooling of 
gas moving the galaxy towards late types (disk dominated). 

2.3 The stellar IMF in quiescent and burst modes of star 
formation 

An important feature of the Baugh et al. model is that stars are as- 
sumed to foiTn with different stellar initial mass functions (IMFs) 
in different modes of star formation. Stars either form quiescently 
in galactic disks or in bursts triggered by some galaxy mergers, as 
explained in the subsection above. In the case of quiescent star for- 
mation, a standard so lar neighbourhood IMF is adopted, following 
the form proposed bv lKennicuttI ( Il983l) . where the number of stars 
produced depends on mass as dN/d\nm oc m~^, with x — 0.4 
for m < IMq and a; = 1.5 for m > IMq. In the case of merger- 
driven starbursts, a top-heavy IMF is assumed where x — 0. The 
use of a top-heavy rather than solar neighbourhood IMF has three 



consequences; i) it increases the amount of radiation emitted in the 
UV, as a larger proportion of massive stars is produced, ii) there is 
higher yield of metals (p = 0.15 compared with p — 0.023 for 
the Kennicutt IMF), leading to more dust and iii) the fraction of 
mass turned into stars which is recycled to the interstellar medium 
is larger (R — 0.91 compared with R = 0.41 for the Kennicutt 
IMF). F or a discuss i on of obser vational support for a top-heavy 
IMF see lLacev et all j2008l,l20I(]|) . 

Finally, in closing this section we remind the read er that an al- 
temative galaxy formation scenario was presented by iBower et al.l 
(2006), who invoke a different mechanism to suppres s the forma- 
tion of galaxies in massive haloes than that used bv iBaugh et al.l 
,Bower et al. also invoke a standard IMF in both modes of star 
formation. Whilst this model gives a good match to the observed 
evolution of the K-band luminosity function and the inferred stel- 
lar mass function, it greatly underpredicts the number of submil- 
limetre galaxies and gives a poorer match to the LBG luminos - 
ity function than the Baugh et al. model (see iLacev et alj|201 ih . 
For furt her discussion of the d i fferences between these two mod- 
els, see 'Almeida et alj ( l2007h . iGonzalez-Perez et al.l ( l2009h and 
iGonzalez et al. ( ,2009i) . 



3 LBG FORMATION HISTORIES AND THE 

PROBABILITY OF HAVING A LBG PROGENITOR 

In this section we present some basic predictions for LBGs in 
the Baugh et al. model. We first compare the model predictions 
for the luminosity function in the UV at z = 3, z = 6 and 
2 = 10 with observations, and use this to set up definitions of 
bright and faint LBGs. (For a more d etailed comparison with ob- 
served LBG luminosity functions see lLacev et alj|20lTI .) We then 
present some illustrative galaxy formation histories for present day 
galaxies which had LBG progenitors. Finally we show the proba- 
bility that a present-day galaxy had an LBG progenitor and that a 
present-day halo hosted an LBG. 

3.1 The UV luminosity function: defining bright and faint 
Lyman-break galaxy samples. 

As outlined in the Introduction, the Lyman-break technique relies 
on colour-colour selection to identify galaxies within a particular 
redshift range, which is set by the redshifted Lyman-break spectral 
feature falling between two of the filters used to image the galaxies. 
In the model we can make predictions for the galaxy population at 
any desired redshift by construction and so there is no need to apply 
any colour selection. Hence, we will use thefarUV(1500A) lumi- 
nosity as our criterion to identify samples of Lyman-break galax- 
ies. We focus our attention on three redshifts, z = 3, z — 6 and 
z = 10, whi ch are representative of the range covered by c urrent 
LBG studies ('Steidel et al." 2003l : lBouwens et alj|2007l . l201 ih . 

[Baugh et al. (2005) presented predictions for the rest-frame 
UV luminosity function at z = 3. Here we revisit this comparison, 
addi ng new observatio nal data at 2 = 3 and going to 2 = 10 (see 
also lLacev et al.ll201 ll) . In Fig.[T]we show the predicted rest-frame 
1500A luminosity function at 2 = 3, 2 = 6 and z — 10, including 
the effects of dust extinction, calculated in a self-consistent way 
using the predicted scale lengths of the disk and bulge components 
and the metallicity of the cold gas. The model matches the obser- 
vational data at 2 = 3, to within the scatt er between datasets, is in 
good agreement with the measurement by Bouwens et alJ (l2007h at 
2 = 6, and is consistent with the tentative measurement and upper 
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Figure 2. Left panels: galaxy merger trees for different z = galaxies. From top to bottom, the stellar masses of the galaxies at z = are M» = 
2.5 X 108,6.6 X 10^" and 2.1 X W^H-'^Mq respectively. The size of the symbol is proportional to the stellar mass as indicated by the key in the top left 
panel. Green circles are drawn for galaxies with Luy > O.lL^y, red circles for Ljjv > L'^^ and blue if Luv < l^ig^t panels: evolution with 

redshift of the galaxy on the main progenitor branch. Top subpanel: evolution of the bulge to total stellar mass ratio (B/T). The dashed vertical line indicates 
the redshift at which the central galaxy (black line) becomes a satellite galaxy (plotted using a green line) after a halo merger. Main light subpanel: the symbols 
show the evolution of the stellar mass (blue), cold gas mass (green), host dark matter halo mass (black) and the stellar mass formed in ongoing bursts (red) in 
units according to the left-hand label. The red line shows the evolutio n of the ratio Ljjy / L*^^ with redshift, in units according to the light-hand axis. Note 
that L'^y is the observationally determined reference value at 2 = 3 iSteidel et all 199^ . 
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Figure 1. The galaxy luminosity function in the rest-frame UV (1500A) at 
2 = 3 (top panel), z = 6 (middle panel) and z = 10 (bottom panel). 
The symbols with error bars show observational estimates taken from 
lAmoutsetai]j2005l) and Reddy & Steidel (2009) at 2 = 3, Bo uwens et all 
)2007l) at 2 = 6 and lfiouwens et alj <2011l) at 2 = 10 as indicated by the 
legends. The predicted luminosity functions include the effect of dust ex- 
tinction and the observational estimates have not been corrected to remove 
extinction. The vertical red lines indicate L^y and O.IL^^- (defined at 
z=3) which we use to construct samples of bright and faint LBGs respec- 
tively. 



limits estimated by iBouwens et alj ilOl ih at z = 10. Following 
the corrmion practice in obsei^ational studies, we define a charac- 
teristic luminosity L^v using the observed position of the break 
in the UV luminosity function at z = 3. We define a galaxy as 
a bright LBG if its UV luminosity is brighter than L^jy and as a 
faint LBG if its UV luminosity exceeds O.lLyy (a sample which 
includes the bright LBG sample, but which is dominated in number 
by galaxies close to the luminosity cut). We apply this definition 
at all redshifts. Note that L^v is taken to be the characteristic lu- 
minosity at 2 ~ 3 found by Steidel et al.i d 19991) (corresponding to 
Mliy = —20.3 + 51og/i in our cosmology). In Fig. [T] we show 
where the O.lL^y and L^y limits lie using vertical lines. 

We note in passing that there is substantial evolution in both 
the observed and predicted UV luminosity functions between z = 
3 and z — 6. The observed characteristic luminosity L^y^^g 
at z — 6 (M*jy , ^fi = —19.5 -I- Slogft in our cosmology, 
iBouwens et alj|2007r) is nearly a magnitude fainter than the z = 3 
value of I/j/y. The number density of galaxies seen at the z — 3 
value of I/f/v drops by around a factor of 5 between z = 3 and 



z = 6. At z = 10, the abundance of galaxies with this luminosity 
is predicted to be several hundred times lower than at 2 = 3. A 
comprehensive study of the evolution of the LBG luminosity func- 
tion, including a more detailed com parison with observational data, 
is presented in lLacev et al.l dZOllh . 



3.2 Example formation histories of model galaxies with 
Lyman-break progenitors 

As examples of the different galaxy formation and merger histories 
which can produce bright LBGs, we plot in the left panels of Fig.|2] 
galaxy merger trees for three present-day galaxies. The mass of 
the galaxies increases down the page. The trees are constructed by 
running GALFORM with many output redshifts. We stop plotting 
the trees at 2 = 6 or when a branch falls below a minimum galaxy 
mass (in stars and cold gas) of W^'h^^ Mq. A branch or trunk of 
the tree is plotted at each output redshift using a circle. The size of 
the circle is proportional to the stellar mass of the galaxy and the 
colour reflects the type of galaxy: green for a faint LBG, red for a 
bright LBG and blue otherwise. 

The galaxy trees are plotted in the following way. Stepping 
back in time from z = 0, at each merger we plot the most mas- 
sive branch on the left and the other branches to the right of this. 
The main progenitor branch is hence the leftmost plotted branch 
traced back from z = in this way. Note that at z > 0, the main 
progenitor branch does not necessarily represent the most massive 
progenitor across the whole tree at a given epoch. The largest pro- 
genitor at any time could jump from one branch of the galaxy tree 
to another, so its selection would not necessarily define a smooth 
path back in redshift. We have chosen to plot examples which have 
at least one LBG in the most massive (leftmost) progenitor branch. 

In the first example (Fig. |2] top left), we show the galaxy 
merger tree for a galaxy which at z — has a stellar mass of 
A/* = 2.5 X lO^/i~^M0 and has a bright LBG progenitor at 
2 = 6. In the second example (Fig. [2] middle left), we show a 
galaxy with M* = 6.6 x W^°h~^MQ at z = with a bright 
LBG progenitor at both z — 6 and z = 3. In the third example 
(Fig. [2] bottom left), we show the galaxy merger tree for a galaxy 
with A/* = 2.1x lO"/i"^M0 atz = with a bright LBG progen- 
itor at z = 6. Note that these trees are purely illustrative examples, 
chosen to show the range of complexity of the trees, and are not 
intended to be a statistically representative sample. 

The right-hand panels of Fig.|2]show the evolution of various 
properties of the galaxy in the most massive progenitor branch. The 
upper inset or subpanel in each case shows the bulge-to-total stel- 
lar mass ratio. If present, a vertical line in this panel marks the 
epoch when the most massive progenitor ceased to be a central 
galaxy and became a satellite. In the examples plotted here, this 
only happens in the case of the least massive galaxy tree shown 
(upper right Fig. |2j, at z ~ 1.1. The main right-hand panel shows 
the UV luminosity of the most massive progenitor, in units of L'^y 
(red line; see right axis for units) and its stellar mass (blue sym- 
bols), cold gas mass (green symbols), stellar mass formed in bursts 
(red symbols) and host dark halo mass (black symbols; see left axis 
for units). In the case of the least massive galaxy (Fig.|2]top right), 
the most massive progenitor experiences a burst at z = 6 which 
makes it a bright LBG. The stellar mass added to the bulge during 
the burst increases the bulge to total stellar mass ratio, and the cold 
gas supply is depleted. Around z ~ 5.2 there is a another burst 
of star formation, triggered by a major merger, which leads to a 
large increase in the bulge-to-total ratio of the galaxy from 0.3 to 
0.9. However, there is less star formation associated with this burst. 
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Figure 3. The triggering mechanism of the LBG phase. For each model 
LBG, we trace its history back in time from the redshift at which it is 
identified, to detemiine whether the UV luminosity is dominated by stars 
formed in a burst triggered by a minor merger (blue), a burst triggered by 
a major merger (red) or stars formed quiescently (green). The top panel 
shows the fractions for these different triggering mechanisms for bright 
LBGs (Ljjv > L^y) identified at z = 3 and 6, while the bottom panel 
shows the fractions for faint LBGs (Ljjv > O.lL'^y ) identified at z = 3, 
6 and 10. (No bright LBGs are predicted at z = 10.) 



due to the depletion of the gas reservoir prior to the burst, and the 
most massive progenitor is only a faint LBG. Another burst hap- 
pens at z ~ 3.5, which makes the most massive progenitor galaxy 
a faint LBG but with relatively little change in the B/T ratio. The 
B /T ratio gradually declines as quiescent star formation adds mass 
to the disk. Finally, after a steady increase in mass, the host halo 
is accreted by a much more massive halo at 2 = 1.1, when the 
most massive progenitor becomes a satellite galaxy. The present- 
day galaxy is a satellite in a halo of mass lO^^h-^MQ. The next 
most massive galaxy (middle panels Fig.|2} displays two epochs at 
which the most massive progenitor is a bright LBG (z ~ 6 and 
z ~ 3.6). The most massive progenitor branch in the most mas- 
sive galaxy (bottom panels of Fig. [2} again shows several instances 
when the progenitor is a bright LBG. It undergoes a major merger 
at z « 1.5 which increases the stellar mass by 30% (barely vis- 
ible in the figure), triggers a bright LBG phase, and changes the 
morphology from B/T ^ to B/T 1. 

It is interesting to note that in all cases shown, the stellar mass 
of the most massive progenitor branch (shown by the blue points in 
the right-hand panels) increases steadily with redshift, without dis- 
playing any sudden large jumps. This implies that the bursts of star 
formation in these examples are not responsible for large changes 
in the stellar mass of the galaxy. This is a result in part of the large 
recycled fraction (R = 0.91) for stars forming with the top-heavy 
IMF, together with the strong feedback from supernova explosions. 



3.3 The triggering and duration of the Lyman-break phase 

We can use the formation histories of model LBGs to answer some 
important questions. The first of these is: what is the mechanism of 
the star formation reponsible for the galaxy being seen as an LBG 
at a particular redshift? In the present model, the possibilities are: a 
starburst triggered by a minor galaxy merger; a starburst triggered 
by a major galaxy merger; and quiescent star formation in a galaxy 
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Figure 4. The distribution of durations of the LBG phase. The top panel 
shows the normahzed distribution of durations for blight LBGs (Luy > 
L^y ) identified at 2 = 3 and 6, while the bottom panel shows this for 
faint LBGs (Luy > O.IL^^) identified at 2; = 3, 6 and 10. The vertical 
aiTows show the median duration for each LBG sample. 



disk. To determine which of these dominates in a particular case, 
we move back in time through the galaxy formation history, start- 
ing at the redshift at which a galaxy is identified as an LBG. If 
at any time during the cuiTent LBG phase the burst and quiescent 
SFRs satisfy the condition SFRquicscont > 4SFRburst then we 
classify the LBG phase as being produced by quiescent star forma- 
tion, otherwise we find the first galaxy merger (going back in time) 
and identify whether it is a minor or major merger. We include the 
factor 4 in the condition on the SFRs because stars forming with 
the top-heavy burst IMF produce roughly 4 times higher UV lu- 
minosities than the same mass of stars formed with the quiescent 
IMF. The results of this classification exercise are shown in Fig. [3] 
for bright LBGs (Luv > L'^y) identified at z = 3 and 6 (top 
panel), and for faint LBGs {Luv > O.lL^v ) identified at 2 = 3, 
6 and 10 (bottom panel). The fractions shown there are calculated 
for a volume-limited sample of LBGs at the identification redshift. 
We see from the figure that LBGs are predicted to be dominated 
by starbursts triggered by minor galaxy mergers (around 60-80% 
of cases) in all of the cases examined here. Minor mergers domi- 
nate over major mergers by a factor 2-A. Quiescent star formation 
is predicted to account for at most a few percent of LBGs, with the 
exception of faint LBGs at 2 = 3, where this fraction increases to 
around 20%. We also see that the quiescent fraction decreases with 
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increasing redshift for both bright and faint LBGs, and is much 
smaller for bright LBGs than for faint LBGs at the same redshift. 
These results about the dominance of bursts over quies cent star for- 
mation in LBGs are in agreement with the analysis bv lLacev et alj 
( l201lh . 

A second question we can answer is: what is the time duration 
of the LBG phase? For bright LBGs, we define this as the time 
during which Luv > ijjy in the current LBG phase, while for 
faint LBGs we define the duration as the time during which Luv > 
O.lLlfY. We show the distribution of LBG durations for volume- 
limited LBG samples in Fig. |4l where the top panel shows results 
for bright LBGs identified at z = 3 and 6, and the lower panel faint 
LBGs identified ai z — 3, 6 and 10. For bright LBGs, the median 
duration is around 20 Myr at both z — 3 and z — 6. For faint 
LBGs, the durations are generally larger, as expected, decreasing 
from around 60 Myr at z = 3 to 30 Myr at 2; = 6 and 20 Myr at 
z — 10. The short durations for bright LBGs reflect the relatively 
short timescales of the bursts responsible for most of them. 

3.4 The fraction of present-day galaxies and hales with an 
LBG progenitor 

Having shown some illustrative examples of galaxy merger trees 
with LBG progenitors in the previous subsection, we now perform 
a statistical study and compute the probability that a present-day 
galaxy had a progenitor that was an LBG or that a present-day halo 
had a progenitor halo which hosted an LBG. 

The predictions for the fraction of present-day galaxies that 
are descendants of LBGs are presented in Fig. |5] with the top pan- 
els showing results for descendants of bright LBGs at 3 < 2 < 4 
and 6 < 2 < 7, and the lower panels results for descendants of faint 
LBGs at 3 < 2 < 4, 6 < 2 < 7 and 10 < 2 < 11. (The model 
does not predict a significant number bright LBGs at z — 10.) In 
the left-hand panels we plot the fraction of present-day galaxies, as 
a function of their stellar mass, which are predicted to have at least 
one bright LBG progenitor in the specified redshift range. In the 
right-hand panels we show the fraction of present-day halos hav- 
ing progenitor halos which hosted an LBG in the specified redshift 
range. We see that the probability to have an LBG progenitor in- 
creases with increasing present-day mass (both stellar and halo), 
but decreases with increasing LBG redshift 

Considering first the bright LBGs, we find that galaxies with 
stellar masses of 6.3 x IO^/i^^Mq, 6.1 x lO^^/i'^M© and 
1.8 X lO"/i~^M0 are predicted to respectively have 10%, 50% 
and 90% probabilities of having at least one bright LBG progenitor 
in the range 3 < 2 < 4. The corresponding stellar masses for bright 
LBGsat6 < 2 < 7 are 4-15 times larger, being 9.0x 10^° /I'^M©, 
2.7 X lO"/i-^M0 and 1.3 x lO^^fc -^Mp respectively. A Mil ky 
Way mass galaxy (4 x 10^°/t~^Mf7^. lBinnev & Tremainell2008l) is 
predicted to have a 40% or 3% probability of having a bright LBG 
progenitor at3<2<4or6<2:<7, respectively. The right 
hand panels of Fig. |5] show that 50% of present-day haloes with 
a mass greater than 4.3 x 10^^ Mq are predicted to have halo 
progenitors that hosted bright LBGs at 3 < 2 < 4. The correspond- 
ing mass for bright LBGs at 6 < z < 7 is an order of magnitude 
larger, 4.8 x W^^H-'^Mq. 

Considering now the faint LBGs, we find that 10%, 50% 
and 90% of galaxies with present-day stellar masses of 1.3 x 
10*/i"^Mo, 2.6 X lO^h-'^MQ and 2.2 x lO^°/i"^M0 are pre- 
dicted to have had faint LBG progenitors at 3 < 2 < 4. These 
masses increase steeply to 1.1 x W^' h^^ AIq, 1.6 x lO^'^h'^ Mq 
and 7.7 x IO^^/i'^Mq for faint LBGs at 6 < 2 < 7, and 



steeply again to 5.3 x W^^h'^MQ, 2.0 x lO^/i'^M© and 
4.2 X 10" ^"^M0 at 10 < 2 < 11. A Milky Way mass galaxy 
is predicted to have a 95% or 70% probability of having a faint 
LBG progenitor at3<2<4or6<2<7 respectively. 
This probability drops to 7% for a 10 < 2 < 11 faint LBG 
progenitor Expressed in terms of halo masses (right-hand panels 
of Fig. [Sj, we find that 50% of halos today with masses larger 
than 2.6 x 10^^h~^MQ are predicted to have halo progenitors 
that hosted faint LBGs at 3 < 2 < 4, while for 6 < 2 < 7 
this increases to 1.4 x 10^^ Mq. Only for masses greater than 
2.7 X W^^h~^M@ are 50% of haloes predicted to have had at least 
one faint LBG progenitor atl0<2<ll. 



4 THE PROPERTIES OF LBGS AND THEIR 
DESCENDANTS 

Having defined samples of bright and faint LBGs, we now com- 
pare the model predictions for their stellar mass, host halo mass 
and morphologies at redshifts 2 = 3, 2 = 6 and 2 = 10. The 
semi-analytical model allows us to follow the subsequent evolution 
of the LBGs, so we also present the predictions for the analogous 
properties of the descendants of the LBGs. We also consider what 
is the dominant mode of star formation - burst or quiescent - in 
LBGs. Note that at z — 10, the model does not predict a significant 
population of bright LBGs (see Fig.[Tll. When we construct the dis- 
tribution of properties of the descendants of LBGs, each descendant 
galaxy is counted once (with a weight depending on its abundance), 
regardless of the number of LBG progenitors it possesses. 

4.1 Stellar masses 

The model predicts that the stellar masses of LBGs identified at 
a given UV luminosity decrease with increasing redshift. On the 
other hand, the stellar masses of their present-day descendants 
show the opposite trend, increasing with increasing redshift. This 
behaviour is seen for both bright and faint LBGs, as we now discuss 
in detail. 

We plot the stellar mass distributions of LBGs and of their 
present-day descendants in the left and right panels of Fig. (6) The 
top panels show results for bright LBGs (Luv > Luy) at 2 = 3 
and 2 = 6, while the lower panels show results for faint LBGs 
(Luv > O.lLuv) at 2 = 3, 2 = 6 and 2 = 10. 

We see that the median stellar mass of bright LBGs at 2 = 3 is 
A/* = 1.3 X 10'-^ H^^Mq, which is 6 times larger than that of their 
faint counterparts at this redshift, and ~ 2.5 times larger than the 
median mass of bright LBGs at 2 = 6. The median stellar mass of 
the 2 = descendants (right panels) of these two LBG populations 
(2 = 3 and 2 = 6) are predicted to be M* = 5 x lO^°/i"^M0 
for the bright LBGs at 2 = 3 and M, = 1 x lO"/i"^A/0 for 
the bright LBGs at 2 = 6. We separate the galaxies into bulge- 
dominated (if the bulge to total stellar mass ratio B/T > 0.5) and 
disk-dominated galaxies (i.e. B/T < 0.5). Bright LBGs at 2 = 3 
and 2 = 6 are predicted to be typically (60%) disk-dominated. 
However, while the descendants of bright LBGs at 2 = 3 are also 
mainly disk-dominated, descendants of bright LBGs at 2 = 6 are 
mainly bulge-dominated. 

The median stellar mass of faint LBGs at 2 — 3 is M* = 
2 X lO^h-^MQ, which is an order of magnitude larger than the 
median mass of faint LBGs at z — 6. Faint LBGs at 2 = 10 are 
25 times less massive than faint LBGs at 2 = 3. The median stellar 
masses of the descendants of the 2 = 3 and 2 = 6 populations are 
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Figure 5. The predicted fraction of present day galaxies (left) and haloes (right) with either bright LBG progenitors (i.e. Luv 
or faint LBG progenitors (i.e. Lyy > O.lL^y ) in bottom set of panels. Left panels: the fraction of present-day galaxies as a function of their stellar mass 
which are predicted to have one or more LBG progenitors. Right panels: the fraction of present-day halos as a function of mass that are predicted to have 
hosted at least one LBG in a progenitor halo. Fractions are shown for LBG progenitors at3<z<4, 6<2<7 and 10 < 2 < 11 (the last for faint 
LBGs only), as labelled. Descendant galaxies are separated into bulge-dominated {B/T > 0.5, red) and disk-dominated {B/T < 0.5, blue). The black line 
corresponds to the case of all LBG descendants and is the sum of the blue and red lines. 
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Figure 6. The stellar mass distribution of LBGs (left) and their z=0 descendants (right). The top set of panels show predictions for bright LBGs {Luv > 
L^^^), and the bottom set of panels faint LBGs (L[/v > O.IL^^^). Results ai'e shown for 2: = 3, z = Gandz = 10 (the last for faint LBGs only) as indicated 
in each panel. The median stellar mass value is given in the label and is indicated by an arrow. Galaxies are separated into bulge-dominated (B/T > 0.5, 
red histogram) and disk-dominated galaxies (B/T < 0.5, blue histogram). The percentage of galaxies in each class is labelled in the same colour as the 
corresponding histogram. The red and blue distributions add to give the total stellar mass distribution (black histogram). 
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predicted to be similar (Af, (3 - 5) x lO^/i^^M©). The median 
descendant mass of a. z = 10 faint LBG is an order of magnitude 
larger than for faint LBGs identified at z = 3 or z = 6. The descen- 
dants of LBGs show a broader distribution of stellar mass compared 
with the LBGs themselves. We also find that faint LBGs at 2: = 6 
and z = IQ are predominantly bulge-dominated systems (75% are 
predicted to have B /T > 0.5). At z = 3 there is a more balanced 
population, with almost half of the faint LBGs predicted to be disk- 
dominated systems. The descendants of faint LBGs identified at 
2 = 3 and 2 = 6 have similar morphological mixes. They are 
predicted to be mainly disk-dominated systems (with 70% having 
B/T < 0.5). The descandants of 2: = 10 faint LBGs are equally 
likely to be disk or bulge dominated. Among the descendants, those 
of lower mass (M < 10^^ Mq) are typically predicted to be 
disk-dominated systems, whereas higher mass LBG descendants 
(M > 10^^ h^^ AIq) are mostly bulge-dominated systems in the 
model. 

Stellar masses of LBGs have been esti mated photometrically 
in a number of ob s ervational studies ( e.g Papovich et al. I I2OOII : 
IShaplev et alj 1200 ll: IVerma et al.l l2007l : Istarket al .1120091) . These 
are all based on fitting model S EDs to multicolour photometry. As 
discussed in iLacev et al.U201lh . these photometric mass estimates 
are fraught with difficulty, since they depend sensitively on the as- 
sumed IMF, and also on assumptions about dust extinction, star 
formation history and metallicity. Since our model galaxies con- 
tain a mixture of stars formed with quiescent and burst IMFs, while 
the observational estimates all assume a solar neighbourhood (typ- 
ically Salpeter) IMF, a direct comparison between stellar masses 
predicted by the model and values estimated observationally from 
photometric data from is not meaningful. Instead. lLacev et aki made 
a direct photometric comparison between the models and obser- 
vations. They plot the observer-frame flux at 3.6/im against rest- 
frame far-UV flux, and find, using the same model as we do here, 
that the predictions come to within a facto r of 3 of matching the ob- 
served relation found by Stark et al. ]([2009) for LBGs at 2 = 4 — 6. 
At the redshifts of the LBG samples, the 3.6/^m flux probes the op- 
tical rest-frame, which is determined by the recent star formation 
history of a galaxy as well as its stellar mass. We will investigate 
this issue in more detail in a future paper. 



4.2 Host halo masses 

The model predicts that the host halo masses of LBGs and of their 
present-day descendants show trends with LBG redshift in the same 
sense as for the stellar masses, i.e. host halo masses for LBGs of 
a given UV luminosity decrease with increasing redshift, while the 
host halo masses of their descendants present-day increase. We now 
discuss this in detail for bright and faint LBGs. 

We plot the predicted distributions of host halos masses of 
LBGs and of their present-day descendants in the left and right 
panels of Fig. |7] The upper panels show results for bright LBGs 
at 2 = 3 and 2 = 6, while the lower panels show results for faint 
LBGs at 2 = 3, 2 = 6 and 2 — 10. 

We see that bright LBGs at 2 = 3 live in halos with a median 
mass of 3 X W^^ h"^ Mq, while bright LBGs at 2 = 6 reside in 
halos with a median mass of 2 x W^^ Mq. The descendants 
of bright LBGs identified at 2 = 3 are hosted by halos with a me- 
dian mass of 1.6 x 10^^h~^ Mq, and the 2 = descendants of 
bright 2 = 6 LBGs are hosted by halos with a median mass of 
5 X 10"/i~^Mq. 

Faint LBGs at 2 = 3 live in halos with a median mass of 
8 X 10^" Mq, which is 4 times larger than the halos which host 



faint LBGs at 2 = 6. The distribution of host halo mass for faint 
LBGs at 2 = 10 is similar to that at 2 = 6. The descendants of 
faint 2 = 3 LBGs are hosted by halos with a median mass of 2.5 x 
10^^ AIq, and the descendants of faint 2 = 6 LBGs are hosted 
by halos with a median mass of 1.3 X 10"/i"^A/Q. Halos which 
host the descendants of faint LBGs have a broader distribution in 
mass than halos which host the progenitors of bright LBGs. 

The host halo masses of LBGs can be observat ionally con- 
straine d through measurements of their clustering. iLacev et al] 
( I2OI ih show predictions for the linear clustering bias b of LBGs 
as a function of luminosity and redshift for the same model as used 
here. For bright LBGs, they find 6 ~ 2 at 2 = 3, increasing to 
& ~ 4 at 2 = 6. As they discuss, these predictions seem broadly 
consistent with observa tion al estimates of the clustering bias by 
lAdelberger elail ( Il998h and lOuchi etHl ( |2004 . 

Earlier calculations with a different semi-analytical model 
suggested that LBGs reside in more massive haloes (» 2 x 
10^^ H-^Mq) at 2 = 3 teaugh et al.lll998l : lGovemato et al.lll998h . 
Ho wever, this earlier m odel had several important differences from 
the lBaugh et alj ^20051) model used in the present paper. It did not 
include the effect of dust extinction on the UV luminosity of the 
model galaxies, it invoked a much stronger supernova feedback 
in low circular velocity haloes, it had much less star formation in 
merger-driven starbursts at high redshift, and it did not assume a 
top-heavy IMF in bursts. The last three of these differences mean 
that in the Baugh et alj feOOSh model, objects can be bright in the 
UV (and t hus be classified as LBGs) while having lower masses 
than in the lBaugh et"ai] ( ll998h model. 



4.3 Morphology 



In §4T|we showed that the model predicts a mixture of morpho- 
logical types (as defined using the bulge to total stellar mass ratio, 
B /T) for LBGs and their descendants. Fig.[8]shows the distribution 
of bulge-to-total stellar mass ratio (B/T) for bright and faint LBGs 
at 2 = 3, 2 = 6 and 2 = 10, and for their present-day descendants. 
We see that faint LBGs are typically more bulge-dominated in stel- 
lar mass than bright LBGs at the same redshift. However, for their 
descendants, this trend is reversed, with the descendants of bright 
LBGs being more bulge-dominated. There is also a tendency for 
LBGs at a given luminosity (whether bright or faint) to become 
more bulge-dominated with increasing redshift, and the same trend 
is also seen for their present-day descendants. 

The observed morphologies of LBGs have mostly bee n stud- 
ied by imaging in their rest-frame UV light. Lacey et alj ( 1201 ih 
present predictions for the B/T ratios of LBGs as seen in rest- 
frame UV luminosity for the same mo del as used here, and com- 
par e with the observa tional studies by iRavindranath et akl ( |2006|) 
and iLotz et aD (|200^. These studies find that for bright LBGs at 
2 ^ 3 — 5, ~ 30% are bulge-like, ~ 10 — 30% a re major mer gers, 
and ~ 40 — 60% are disk-like or minor mergers. iLacev et akl con- 
clude that the model predictions are qualitatively consistent with 
these observational results. 



4.4 Mode of star formation in LBGs 

As discussed in Section 2, the model assumes two modes of star 
formation, quiescent star formation in disks and star formation in 
bursts triggered by galaxy mergers. Quiescent star formation pro- 
duces stars with a standard solar neighbourhood IMF, whereas in 
burst star formation, a top heavy IMF is adopted. This difference 
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Figure 7. Host halo mass distribution of LBGs (left) and tlie liost tialo mass distribution of their z = descendants (right). The top set of panels show bright 
LBGs (LuY > L^y) and the bottom set of panels show faint LBGs (Ljjy > O.lL^y). Predictions are shown for 2 = 3, z = 6 and z = 10 (the last for 
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in the choice of IMF between the two modes of star formation has 
an impact on the abundance and properties of LBGs. In this sec- 
tion we study correlations between the star formation rate, stellar 
mass, halo mass and UV luminosity of galaxies at z — 3, z = 6 
and z = 10. Our goal is to understand the role of the different star 
foimation channels in the production of LBGs and in determining 
their nature. 

In Fig. 121 we plot all galaxies with stellar masses M* > 
10^ Mq at 2 = 3, 6 and 10, showing the relations between 
stellar mass (plotted on the x-axis) and various properties on the y- 
axis: host halo mass (left panel), star formation rate (middle panel) 
and far-UV (1500A) luminosity in units of L^jy (right panel). (We 
choose a lower mass limit of A/, = 10^h~^ Mq because the 
masses of faint LBGs at z = 10 extend down to roughly this mass.) 
The shading in Fig.|9]reflects the square root of the space density of 
galaxies, becoming darker where more galaxies are found. We gain 
further insight by dividing the galaxies according to which mech- 
anism of star formation is dominant (but without the factor 4 used 
earlier in |J33) and whether the galaxies are central or satellites in 
their host dark halo. Contours of a given colour enclose 95% of the 
galaxies passing one of these selection criteria. The red contours 
delimit the region occupied by 95% of the galaxies in which star- 
bursts are the dominant mode of star formation and blue contours 
enclose the quiescent population. Similarly, orange contours show 



where the majority of central galaxies lie and green contours show 
where the satellite galaxies fall. 

Fig. [9] shows that there is a correlation between galaxy stellar 
mass and host halo mass which is driven by central galaxies. Satel- 
lite galaxies in the model, which are mostly quiescent, are respon- 
sible for the large scatter in this relation (as suggested by the green 
contour). There is also a correlation between stellar mass and star 
formation rate. In this case the correlation is largely due to galaxies 
where quiescent star formation dominates (as suggested by the blue 
contours), with the scatter produced mainly by galaxies which are 
undergoing bursts of star formation (see also lLagos et al.lllOlOl ). 

In Fig.[TO]we plot similar correlations, but in this case for faint 
LBGs (L > O.lLlfy ) only. Again, we differentiate between galax- 
ies according to the dominant mode of star formation, and make 
this distinction for the faint and bright LBG samples. Fig.llOlshows 
that most LBGs are produced by b ursts of star format ion, as dis- 
cussed earlier in i]3.3l and also by iLacev et alj i201lh . We again 
see that the fraction of LBGs dominated by burst star formation in- 
creases with redshift. The correlations of halo mass and SFR with 
stellar mass are seen to be better defined, with less scatter, for LBGs 
than for samples containing all galaxies. However, the correlation 
of SFR with stellar mass is also seen to be different (shifted to 
higher SFRs) for LBGs from that for all galaxies. This is because 
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quiescent star fonnation rate and the blue contour shows the converse. The orange contour marks the locus within which 95% of central galaxies are to be 
found and the green contour shows the corresponding distribution of satellites. The labels of the same colour as the contours give the percentage of galaxies 
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the LBG samples are dominated by bursts, which have larger SFRs 
than quiescent galaxies of the same stellar mass. 

We have also investigated what fraction of the stellar mass of 
model LBGs has formed quiescently or in bursts (allowing for the 
recycling of mass from dying stars). We find that LBGs at 2 = 3 
typically formed most of their stellar mass by quiescent star for- 
mation, with 57% of faint LBGs and 73% of bright LBGs having 
formed more than half of their stellar mass quiescently. However, 
at higher redshifts, the reverse is true: at z = 6, 93% of faint LBGs 
and 89% of bright LBGs have built at least half their stellar mass 
through bursts of star formation after merging with other galaxies, 
while at z = 10, 99% of faint LBGs have built at le ast half their 
stella r mass through bursts. (See also the discussion in lLacev et al.l 
I20TI) 



5 WHAT IS THE CONNECTION BETWEEN LBGS AND 
SMGS? 

Samples of galaxies at high redshift have been assembled using 
the Lyman-break technique and through their emission at sub- 
millimetre wavelengths. LBGs are selected on the basis of their 
emission in the rest-frame far-UV. SMGs are detected through their 
rest-frame far-IR/sub-mm emission, which is produced by the rera- 
diation by dust of light absorbed predominantly in the UV. It is 
natural to ask what is the overlap, if any, between these two popu- 
lations at high redshift. A first indication of the answer is suggested 
by the different redshift ranges of the populations. LBGs selected 
using ground-based surveys are restricted to have z > 3. Bright 
SMGs (with 850/^m flux Sgso > 5 mjy), on the other hand, have 
a median redshift of z ~ 2 in flux-limited samples. Hence, to study 
the overlap between the two samples we will be looking at the UV 
luminosities of galaxies in the high redshift tail of the SMG redshift 
distribution. Almeida et al. (2010) performed a similar comparison 
at 2 = 2, and found that ~ 1% of galaxies bright in the UV were 
predicted to also be bright SMGs. 

The multiwavelength nature of our semi-analytical model and 
its self-consistent treatment of dust make it ideally placed to ad- 
dress the overlap between the SMG and LBG populations. Here we 
examine the predicted relation between 850/im flux and UV lumi- 
nosity at z = 3 and z — 6. Following the companion to this paper 
jGonzalez et al.ll201lh . in which we study the evolution of SMGs 
in the model, we define a galaxy as a bright SMG when its 850/im 
flux exceeds 5 mJy and as a faint SMG when its 850/im flux ex- 
ceeds 1 mJy. 

The top panels of Fig. [TT] show the relation between the 
850 fivci flux and the UV luminosity for galaxies at 2 = 3 and 
2 = 6, respectively. There is scatter in this relation, particularly at 
fainter magnitudes. Nevertheless an increase in the UV luminosity 
typically corresponds to an increase in the submillimetre flux. 

The middle panels of Fig. [TT] show the number density of 
galaxies as a function of their 850/im flux, with separate curves 
for faint and bright LBGs at 2 = 3 and 2 = 6, respectively. At both 
redshifts, the median fluxes at 850/im are predicted to be l.OmJy 
for bright LBGs, equal to the faint SMG flux limit adopted here, but 
well below it (~ 0.1 - 0.2 mJy) for faint LBGs. We find that only 
1.4% of the bright LBGs {Luv > L\jy) at 2 = 3 are also bright 
SMGs (5860 > 5mJy). This figure decreases to 0.1% for 2 = 6. 

In the bottom panels of Fig. [TT] we show the predicted num- 
ber density of galaxies as a function of UV absolute magnitude at 
2; = 3 and 2 = 6 respectively, distinguishing between faint and 
bright SMGs. At 2; = 3, the median UV magnitude is quite sim- 



ilar for faint and bright SMGs (~ 0.51/*). At 2 = 6 the median 
UV absolute magnitude for the faint SMGs is predicted to be very 
similar to its value at z = 3. Our model predicts too few bright 
SMGs at z — 6 for us to reliably calculate the median UV abso- 
lute magnitude (as suggested by the redshift distributions plotted in 
iBaugh et al.l2005l) . At 2 = 3, nearly half (42%) of the bright SMGs 
are predicted to also be bright LBGs. 

The model predicts that LBGs at 2: = 10 should be extremely 
faint at 850 /im, with faint LBGs having median fluxes of only 
0.03 mJy. We predict a negligible number of bright or faint SMGs 
at 2 = 10. 



6 SUMMARY 

We have investigated the nature, formation histories and present- 
day descendants of Lyman-break gal axies (LBGs) in the ACDM 
cosmology that are predicted by the iBaugh et al.l ( |2005|) version 
of the GAL FORM semi-analytical galaxy formation model. The 
'Baugh et al.' model incorporates a top-heavy IMF in starbursts trig- 
gered by galaxy mergers, which was found to be necessary to 
explain the number counts and redshift distribution of sub-mm 
galaxies. This model has previously been shown to match the ob- 
served far-UV l uminosity functio n of LBGs over the whole range 
z — 3—10 ( Lacev et alj|201ll) . In this paper we investigate a 
range of other predictions for Lyman-break galaxies from the same 
galaxy formation model. 

We began by investigating the formation histories of LBGs. 
We found that the majority of LBGs are predicted to be starbursts 
triggered by minor mergers of galaxies. Using the observed char- 
acteristic far-UV luminosity at 2: ~ 3, L^jy, to define samples of 
model galaxies as bright (I/(7v > L^y) or faint (Luv > O.lL^v) 
LBGs, we found that > 99% of bright LBGs at 2 = 3 - 6 are pre- 
dicted to be starbursts triggered by minor or major mergers, while 
this fraction is somewhat lower for faint LBGs (but still > 80% 
for 2 = 3 — 10). The importance of starbursts in producing galax- 
ies visible as LBGs is because the far-UV light is dominated by 
young stars, and in starbursts the star formation rate is increased 
over the quiescent value. However, in our model there is also sec- 
ond important effect, that the top-heavy IMF in bursts (as against 
a solar neighbourhood IMF assumed for quiescent star formation) 
further enhances their far-UV luminosities over quiescent galaxies. 
We then examined the duration of the LBG phase in model galax- 
ies, and found that this is quite short, around 20 Myr for bright 
LBGs, and somewhat longer, around 20 — 60 Myr, for faint LBGs. 
The short timescales are a result of the LBGs mainly being star- 
bursts. 

We then computed the fraction of present-day galaxies that are 
predicted to have LBG progenitors. We find that 95% of Milky Way 
mass galaxies had at least one faint LBG progenitor in the redshift 
range 3 < 2 < 4, and 70% at 6 < 2 < 7, but this falls to 7% 
for 10 < 2 < 11. For bright LBG progenitors these fractions drop 
to 40% for LBGs at 3 < 2 < 4 and only 3% at 6 < 2 < 7. 
Performing a similar analysis for dark matter halos, we find that 
50% of present-day halos with masses exceeding 4 x W^^ Mq 
or 5 X W^'^h'^ Mq had progenitor halos hosting bright LBGs at 
3<2<4or6<2<7 respectively. 

We next examined the stellar and halo masses of LBGs, to- 
gether with their morphologies. The model predicts that bright 
LBGs at 2 = 3 have median stellar mass ~ 1x10^ Mq , falling 
to ~ 5 X W^h~^ Mq at 2 = 6, and are typically disk-dominated 
in stellar mass, consistent with their origin in starbursts triggered 
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Figure 11. The relation between 850 fj,m flux and rest-frame far-UV absolute magnitude for galaxies at z = 3 (left panels) and z = 6 (right panels). Top 
panels: the median 850 fj.m flux as a function of the rest-frame far-UV absolute magnitude for galaxies at z=3. The bars indicate the 10 to 90 percentile range of 
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by minor galaxy mergers. As already noted by iLacev et al.l (1201 ih . 
these stellar masses are lower than the values typically inferred 
from observational data by fitting stellar population models to 
broad-band photometry, but the observational estimates rely on as- 
suming a solar neighbourhood IMF, while in our model the light 
from LBGs is dominated by stars formed with a top-heavy IMF 
during a starburst. The faint LBGs are predicted to have stellar 
masses roughly an order of magnitude lower than bright LBGs at 
the same redshift, falling to only ~ 1 x 10^ Mq at z = 10, and 
to be more bulge-dominated than the bright LBGs, especially at the 
higher redshifts. The halos hosting bright LBGs at 2; = 3 and z — 6 
are predicted to have median masses ^ 10";i~^Mq, with roughly 
10 times lower masses for faint LBGs at the same redshifts. 

Following on from this, we then investigated the stellar and 
halo masses and morphologies of the present-day descendants of 
LBGs. The ability to relate different types of high-redshift galax- 
ies to their descendants is one of the unique capabilities of the- 
oretical models of galaxy formation such as ours. We find that 
the descendants of bright LBGs at 2 = 3 have median stellar 
masses ~ 5 x 10^"/i~^Mq (i.e. comparable to the Milky Way), 
and consist of roughly equal numbers of disk- and bulge-dominated 



systems. They typically inhabit galaxy groups, with halo masses 
~ lQ^^ft"^M0. Bright LBGs at 2 = 6 have descendants which 
are somewhat more massive 1 X lO"/i"^M0) and inhabit 
somewhat more massive halos, but are significantly more bulge- 
dominated (~ 75%) on average. Faint LBGs are predicted to have 
present-day descendants which are an order of magnitude lower in 
stellar and halo mass than the bright LBGs at the same redshift, and 
contain a larger proportion of disk-dominated systems. 

We next looked at the relationship between halo mass, stel- 
lar mass, star formation rate (SFR) and far-UV luminosity in the 
model, for all galaxies and for LBGs. We find that LBGs have a 
tighter relationship between halo mass and stellar mass than galax- 
ies as a whole, due to them them being mainly central galaxies. We 
also find a correlation between SFR and stellar mass that is shifted 
to higher SFRs compared to that for all galaxies. This results from 
LBGs being mainly starbursts, which have higher SFRs than qui- 
escent galaxies with the same stellar mass. We also find that al- 
though most LBGs in the model are starbursts, their stellar masses 
are not necessarily dominated by stars formed in bursts. This re- 
sults in part from the large recycled fraction for stars forming with 
the top-heavy burst IMF. LBGs at 2 = 3 (both faint and bright) are 
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predicted to have formed most of their stellar mass quiescently be- 
fore the cuiTent burst. However, this changes with redshift, so that 
LBGs at 2 = 6 have typically formed most of their stellar mass in 
bursts (although not necessarily the current burst). 

Finally we used the multi-wavelength capabilities of our 
model to examine the overlap between LBGs and sub-mm galax- 
ies (SMGs) detected at 850/im, extendi ng to z = 3, 6 and 10 
the comparison previously carried out by lAlmeida et al.l fcoioh at 
2 = 2. Analogously to the LBGs, we define bright and faint SMGs 
as having 850/im fluxes exceeding 5 mjy and 1 mJy respectively. 
At z ^ 3 we are only probing the tail of the redshift distribution of 
bright SMGs, which peaks at 2 ~ 2. We find that only about 1% 
of bright LBGs at 2: = 3 are predicted to be also bright SMGs, and 
this fraction falls to only 0.1% for bright LBGs at 2 = 6, with the 
median 850/im flux being about 1 mJy for both redshifts. This ex- 
plains the difficu lty in observationally detecting the sub-mm emis- 
sion from LBGs JChapman et alj2OO0l) . 
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